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ABSTRACT 

We present the discovery of compact, obscured star formation in galaxies at z ~ 0.6 that exhibit 
> 1000 km s _1 outflows. Using optical morphologies from the Hubble Space Telescope and infrared 
photometry from the Wide-field Infrared Survey Explorer, we estimate star formation rate (SFR) 
surface densities that approach Ssfr ~ 3000 M© yr^ 1 kpc~ 2 , comparable to the Eddington limit from 
radiation pressure on dust grains. We argue that feedback associated with a compact starburst in the 
form of radiation pressure from massive stars and ram pressure from supernovae and stellar winds is 
sufficient to produce the high- velocity outflows we observe, without the need to invoke feedback from 
an active galactic nucleus. 

Subject headings: galaxies: evolution — galaxies: kinematics and dynamics — galaxies: ISM - 
galaxies: starburst 



1. INTRODUCTION 

The central regions of elliptic al galaxies are thought 
to form in compact starbursts (jKormendv et al.l 120091 : 
iHopkins et al.l I2009D . Feedback associated with such 
starbursts can produce outflo ws driven by thermal energ y 
from superno va explosions ([Chevalier fe Cleed [1985), 
stellar winds ([Leitherer et al.lfi992D . and momentum in- 
put from both supern ova ram pressure a nd radiation 
pressure on dust grains ()Murrav et al .1120051) . It has been 
argued that such feedback imposes a limit on the max- 
imum star-for mation rate (SFR) surface density (Ssfr.) 
for starbursts (ILehnert fc Heckmanl 119961 jMeurer et al.l 
[19971: IMurrav et al.ll2005t iThompson et al.l l2005) and the 
maximum stellar surface density for elliptical galaxies 
and star clusters ()Hopkins et al.ll2~010D . 

Galactic winds are ubiquitous in star-forming galaxies 
at all redshifts and generally exhibit outflow velocities 
in the 100-500 km s _1 range, which ca n be attributed 
to th e stellar processes des c ribed above (IHeckman et al. | 
2000] IShaplev et all [2001 [Martini [20051: IRupke et al.l 
20051: IWeiner et al.l I2009L IRubin et al.l 12010ft . Outflows 
with significantly hi gher velocities (l y| > 1000 km s _1 ) 
were discovered by iTremonti et al.l (|2007l ) in a sample 
of massive (M* ss 10 11 M ) post-starburst galaxies at 
z ~ 0.6, and it was suggested that a more energetic 
source suc h as feedback from an accreting superma ssive 
black hole (|Silk fc Rccs 1998: lDi Matteo et al.ft2005j) may 
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be responsible for launching the winds (see iFab 
for a recent review) 
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However, it is also plausible that feedback from a com- 
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pact starburst could expel gas with such large velocities. 
Indeed, there is evidence for a positive correlati on be- 
tween outflow velocity and starburst luminosity ([Martini 
[2005t IRupke et al.l l2005t ITremonti et"aT1 [200%. albeit 
with significant scatter. Furthermore. IHeckman et al.l 
(|2011[) recently found outflows with maximum velocities 
reaching 1500 km s _1 in a sample of local starbursts with 
compact nuclei, and argued that such velocities could be 
explained by a wind launched from r$ ~ 100 pc and 
driven by feedback from massive stars and supernovae. 

In this Letter, we measure sizes and SFRs for a sam- 
ple of massive galaxies at z ~ 0.6 that exhibit > 
1000 km s -1 outflows , expanding on the initial study by 
ITremonti et al.l ()2007f ). We seek to test whether the en- 
ergetic outflows in these galaxies could have been driven 
by feedback from starbursts with very large SFR surface 
densities. Our analysis combines galaxy sizes measured 
with the Hubble Space Telescope (HST) with SFRs and 
stellar masses estimated from Wide-field Infrared Survey 
Explorer (WISE), Spitzer Space Telescope, Sloan Digi- 
tal Sky Survey (SDSS), and Galaxy Evolution Explorer 
(GALEX) photometry. 

2. ANALYSIS 
2.1. Morphologies and Sizes 

We observed 29 galaxies with HST (programs 12019 
and 12272) selected from a parent sample of ~ 10 3 
galaxies at 0.35 < z < 1.0 with post-starburst spec- 
tral features: B or A-star dominated stellar continua 
and moderately weak nebular emission (EW([0 n]) < 
20 A; see C. Tremonti, et al., in preparation for more 
details). The galaxies targeted with HST were those 
with the youngest derived post-burst ages, % urst < 
300 Myr. Therefore this sample has bluer U ~ B colors 
and stro nger emission lin es than typical post-starburst 
samples (|Coil et al.ll2011l) . Our subsequent UV-IR SED 
analysis (see Section I2.2|) reveals significant ongoing star 
formation (SFR > 50 M Q yr" 1 ) in the 14/29 galaxies 
with WISE 22 (im detections, calling into question the 
post-starburst nature of roughly half of the HST sample. 
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Fig. 1.— Left: Observed UV-IR SEDs (A rest = 0.1-15 /^m) 
for three galaxies with the largest SFR surface densities (Egpj^ fa 
3000 Mq yr _1 kpc -2 ). The top two SEDs are offset by 5 and 
10 magnitudes, respectively. We show stellar population fits to 
the A res t = 0.1-3 fira emission (black solid line) and three tem- 
plates for dust emission (M82 starburst, blue dashed line; Arp 220 
star burst, green dot -dashed line; obscured quasar, red dotted line; 
IPolletta ct al. 2007) scaled to match the A ^ s = 4.6 fim band. 
The starburst templates provide reasonable fits to the A Q b s = 12 
and 22 /an WISE photometry, while the quasar template does not. 
Right: HST/WFC3 F814W images (probing A rcst 5000 A) show- 
ing that these galaxies are dominated by a compact nucleus. 



Using the F814W filter on the WFC3/UVIS channel, 
which has 0.04" pixels, we obtained 4 x 10 min exposures 
in a single orbit for each galaxy. The dithered images 
were processed with MultiDrizzl^H to produce science 
mos aics with 0.02" pix els. For each galaxy, we use GAL- 
FIT (jPeng et al.ll2002l ) to model the two-dimensional sur- 
face brightness profile with a single Sersic component 
(characterized by Sersic index n and effective radius r e ), 
using stars in the images to construct the model point- 
spread function (PSF) . In cases where the best-fit model 
returns n > 4, we also fit an n — 4 de Vaucouleurs model, 
yielding a larger r e value (due to the covariance between 
n and r e ); we use these larger effective radii in our anal- 
ysis. 

In this Letter, we highlight the galaxies with the small- 
est r e and largest Ssfr values because such extreme 
starbursts could conceivably produce the high-velocity 
outflows we observe (see Section [3]). We show HST im- 
ages for the three highest Esfr galaxies in Figure [TJ In 
all three cases, the single-component GALFIT model ac- 
counts for > 85% of the total flux. The residuals show 
diffuse emission that is consistent with these systems be- 
ing late-stage galaxy mergers, although we defer a de- 
tailed study of the merger stage to future work. 

For the most compact galaxy (J0905+5759, r e = 
0.013" or 100 pc), we also show the observed one- 
dimensional surface brightness profile in Figure [5] We 
compare to the profiles of six stars in the same image, 
the best-fit de Vaucouleurs model, and a de Vaucouleurs 
model with r e = 0.04" (the native WFC3/UVIS pixel 

7 http:/ /stsdas. stsci.edu/multidrizzle/ 
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Fig. 2. — One-dimensional surface brightness profile for 
J0905+5759, which has the smallest effective radius in the sample. 
The observed profile is shown as the solid black line, and the pro- 
files of six stars in the same image are shown in blue, normalized to 
the same central surface brightness. The best-fit de Vaucouleurs 
profile with r e = 100 pc is shown as a dashed red line, and for 
comparison a broader profile with r e = 0.04" = 290 pc is shown as 
the dotted green line. This galaxy is quite compact, but is more 
extended than a point source. 



size, which corresponds to a physical scale of 290 pc). 
This comparison illustrates that this galaxy, while only 
marginally resolved with an r e that is ~ 20% of the image 
FWHM, is clearly more extended than a point source. 

For such a compact source, there is uncertainty in our 
r e measurement given uncertainties in the model PSF. 
To quantify this, we used TinyTim to generate a model 
PSF that is narrower than the stars in the image (con- 
volving with a FWHM = 0.04" Gaussian, whereas the 
image FWHM is 0.074") and found that this increased 
the r e in the GALFIT model by a factor of two. We also 
fit a two-component PSF+Sersic model, but found that 
the Sersic component dominates the fit, yielding a sim- 
ilar r e . Furthermore, the spectrum of the galaxy shows 
no evidence for an AGN contribution to the optical con- 
tinuum (see Figure [3]) , so there is no clear motivation 
for including an unresolved, point-source component in 
the model. We conclude that our r e estimate is accurate 
within a factor of two. 

2.2. Star Formation Rates and Stellar Masses 

We gathered photometry from the WISE All-Sky Re- 
lease, the SDSS Seventh Data Release, and GALEX Gen- 
eral Release 6. We also obtained 5 x 30 sec dithered 
exposures at 3.6 /im and 4.5 /im for all sources with 
the Infrared Array Camera on the Spitzer Space Tele- 
scope as part of GO program 60145. We used the post- 
basic calibrated data to perform aperture photometry 
on all sources and point-source photometry on sources 
in crowded fields. We show spectral energy distributions 
(SEDs) for the three highest Ssfr galaxies in Figure [TJ 

We estimate IR-based SFRs for the 25/29 galax- 
ies with WISE 12 or 22 /im detections by fitting 
iCharv fc Elbazl (|2001[) templates to their 12 and 22 yttm 
fluxes. For the 14/25 galaxies with 22 /im detections, 
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Fig. 3. — Spectra covering A rcst = 2500-5200 A for the three 
galaxies shown in Figure [T] For clarity the top two spectra are 
offset by +10 and +20 units. These spectra are dominated by 
the light of a young stellar population but have relatively weak 
nebular emission lines ([O II] A3727, H/3 A4861, [O III] A5007) and 
strong Mg II AA2796, 2803 absorption arising from the interstellar 
medium. The absence of broad Mg II or H/3 emission lines rules 
out a significant contribution to the optical continuum light from 
an AGN. The panels on the right highlight the region around the 
Mg II doublet in both wavelength and velocity space to illustrate 
the outflow kinematics. The spectrum in the bottom panel has 
sufficient spectral resolution (FWHM 8 km s _1 ) to resolve the 
intrinsic shape of the absorption-line profile, revealing that the gas 
near the centroid velocity (v = —2470 km s —1 ) covers the entire 
galaxy. 



this yields SFRs that agree with those obtained from 
the iRuiopakar n et all (|2012l ) method based on 24 /zm 
luminosity with a scatter of 0.05 dex. Several authors 
have shown that the shape of the IR SEP for star- 
forming galaxies depen ds on Esfr (iRuiopakar n et al.l 
1201 It lElbaz et al.l 120111) , with more compact starbursts 
having larger total-IR (8-1000 /Ltm) to mid-IR (8-24 fim) 
ratios, charact eristic of the most luminous galaxies in the 
local universe (|Rieke et al.M2009h . If we used the most lu- 
minous local templates for the 8/25 sources with SFRs 
in the ULIRG regime (SFR m > 100 M yr -1 ), we would 
obtain SFRs that are larger by 0.5 dex than the values 
we adopt for this Letter. 

We also estimate SFRs and stellar masses based on 
stellar populatio n fits to the A rss t = 0-1- 3 /an SEDs us- 
ing the method of lMoustakas et al.l (|2011| ). For the 14/25 
galaxies with SFRtr > 50 M© yr -1 , there is agreement 
between these UV-based SFR estimates and SFRtr with 
a scatter of 0.32 dex. For an SMC dust law, we find a me- 
dian attenuation of Ay = 0.4 mag. The observed H/3 lu- 
minosities, uncorrected for dust extinction, are typically 
factors of 10-20 x smaller than expected from the UV and 
IR SFRs. This can be reconciled by either strong differ- 
ential dust attenuation (i.e., Ay w 2-3 mag for the line- 
emitting region), escaping ionizing photons from matter- 
bounded H II regions, or a recently quenched starburst 
(t > 5 Myr) with a small ratio of ionizing (A < 912 A) 
to non-ionizing UV photons. 

2.3. Outflow Kinematics and Covering Factors 



We present A rost = 2500-5200 A spectroscopy for 
three high-EsFR sources in Figure [3] based on data 
from MMT/Blue Channel and SDSS (J1506+5402), 
Magellan/MagE (J1341-0321), and Keck/LRIS and 
Keck/HIRES (J0905+5759). The spectra are dominated 
by light from a young (t < 50 Myr) stellar population. 
We highlight the interstellar medium Mg n AA2796, 2803 
absorption lines, which are used to measure outflow 
velocities. At low spectral resolution (e.g., the top 
right panel of Figure [3]) it is not possible to deter- 
mine the intrinsic shape of the absorption line pro- 
file and therefore the covering factor of the outflowing 
gas. However, the Keck/HIRES spectrum of J0905+5759 
(FWHM ss 8 km s" 1 ) reveals that the gas covers the en- 
tire continuum source near the velocity centroid (v = 
—2470 km s -1 ) indicating a galaxy-wide outflow. 

3. DISCUSSION 

The compact sizes (r e w 100 pc) and large SFRs 
(SFR k, 200 Mg) for the three galaxies described above 
imply extremely large SFR surface densities (Esfr ~ 
3000 M yr" 1 kpc~ 2 ). To place these galaxies in con- 
text, we plot Esfr versus stellar mass for the 25/29 
galaxies detected by WISE in Figure 21 We include com- 
parison sam ples of ~ 10 5 star-fo rming galaxies at 0.5 < 
z < 1.5 from lWuyts et al.1 1)201 lh a nd gas-rich me r gers a t 
z < 0.3 including 32 ULIRGs from lVeilleux et al.l (|200l , 
six L yman break analogs w ith dominant central objects 
from lOverzier et al.l (120091) . and the local compact star- 
burst Arp 220 (jKennicuttj 119981 iRodriguez Zaurin et al.1 
l2008h . We also mark the empirical threshold for launch- 
ing w inds (Esfr ~ 0.1 M yr _1 kpc -2 , iHeckmanl 
2002), the 90th percentile limit for the surface bright- 
ness of starbursts over a wide range in redshift mea- 
sured using UV, Ha, far-IR, and radio continuum emis- 
sion (Esfr ~ 25 M yr -1 kpc" 2 for a Chabrier IMF, 
iMeurer et al.l I1997D . and the theoretical limit for a 
starburst limited by feedback fro m radiation pressure 
(Esfr ~ 3000 M^ yr -1 kpc" 2 . iMurrav et all [20051: 
iThompson et al.l 120051 : iHopkins et al.l l2010f ). The most 
luminous, compact starbursts in our sample exhibit SFR 
surface densities that reach the Eddington limit, suggest- 
ing that their growth is being regulated by momentum 
input from massive stars. 

3.1. Constraints on Ongoing AGN Activity 

While the SEDs (Figure [1]) and optical spectra (Fig- 
ure [3]) for our sample indicate that their A ros t = 0.1- 
15 /zm emission is dominated by star formation, it is 
worthwhile to consider the level of ongoing AGN activ- 
ity and its effect on our measurements. We observed 
12/29 galaxies with the Chandra X-ray Observatory 
(proposal ID 11700896, see P. Sell et al., in preparation 
for more details), including the 3/29 galaxies with signifi- 
can t [Ne v] A3426 de tections (an indication of AGN activ- 
ity, iGmTeLinioIii and the most luminous [O ill] A5007 
sources in the sample. Of these 12 galaxies, only two 
were detected with > 4 X-ray counts, and upper lim- 
its for the remainder imply Lx < 10 43 erg s _1 . The 
two significant detections were for galaxies hosting op- 
tical broad-line AGNs; neither was detected at 22 /im 
by WISE or shows evidence for an outflow, and both 
have Esfr < 10 M yr _1 kpc -2 . For the most lumi- 
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Fig. 4. — SFR surface densities and stellar masses for the HST-WISE sample described in this paper (black circles, symbol size pro- 
portional to o utflow velocity), along with samples o f z < 0.3 gas-rich mergers (orange squares, including ULIRGs, Lyman break analogs, 
and Arp 220; IKennicutt 1998; Vcillc uTeFaTI [2003; IRodriguez Zaurin et af] [20081; lUverzier etUI [2009) and z ~ 1 star-forming galax- 
ies (shown with 68%, 95%, and 99.7% contours; IWuvtsetlrtTl2*o"TTrK We mark the empirical threshold for launching winds (dotted 
line, SgpR ~ 0-1 Mq yr _1 kpc~ 2 ; Hcckman 2002), the 90th-percentile starburst intensity limit from Mcurcr et al. (1997) (dashed line, 
Ssfr ~ 25 Mp) yr —1 kpc~ 2 ), and the Eddington limit from ra diation pressure on dust grains (solid line, SgpR 3000 Mq yr _1 kpc - 2 ; 
Murray et al. 2005; Thompson ct al. 2005; Hopkins ct al. 2010). The representative error bar in the top-right portion of the plot corre- 
sponds to uncertainties of 0.3 dex in SgpR an d 0.2 dex in stellar mass. Our HST— WISE sample overlaps with the region characterized by 
gas-rich mergers, and extends to very large SFR surface densities near the Eddington limit, suggesting growth that is limited by momentum 
injection from massive stars. 



nous [O in] source, J1506+5402, the observed [O iii]/X- 
ray ratio implies X-r ay absorption by a factor of ~ 10 
(jHeckman et al.f2005l ). However, even if the X-ray atten- 
uation were a factor of 100, typical for l ocal Compton- 
thick AGNs (jDiam ond-Stan ic et al.ll2009fl. the e xpected 
mid-IR AGN contribution (jGandhi et all I2009D would 
still be <30% of the observed mid-IR luminosity. Only 
1/29 galaxies has WISE mid -IR colors that w ould iden- 
tify it as an AGN candidate ()Stern et al.ll2012t ). and this 
particular galaxy has Esfr < 50 M yr" 1 kpc -2 . We 
conclude that the bolometric output of the galaxies in 
our sample is dominated by star formation and that our 
results regarding large Ssfr values are not affected by 
AGN contamination. 



3.2. The Outflow Launching Mechanism 

Are the high-velocity outflows observed i n these galax- 
ies pr oduced by a compact starburst? iMurrav et al.1 
(|2011[) argued that massive star clusters with large 
gas surface densities are the ideal launching point for 
galactic-scale outflows driven by radiation pressure, and 
that the outflow velocity should scale with escape veloc- 
ity of the most massive star clusters in a galaxy. For 
our sample, if one assumes that the spatial extent of 
the stellar mass is comparable to that of the rest-frame 
y-band light (see Section T3.3p . then half of the stellar 
mass (M* ~ 10 11 M©) will be within the effective radius 
(r e ~ 100 pc). Such a compact stellar population would 
have an escape velocity comparable to the > 1000 km s _1 
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This argument, combined with the fact that we observe 
galaxies with significant dust-obscured star formation 
and Ssfr values near the Eddington limit, suggests that 
momentum input from massive stars in the form of radi- 
ation pressure is a viable mechanism for launching these 
outflows. 

In addition to radiation pressure, we also expect signif- 
icant momentum flux from stellar winds and supernovae. 
For example, a starburst with SFR « 200 M Q yr" 1 would 
have radiation pressure Ltoi/c «3x f0 35 dyne and ram 
pressure p « 5-10 x 10 35 dyne from ste l lar winds and 
supernovae dLeitherer et al.l 119921 Il999t iVeilleux et al.l 
120051 ). Furthermore. iHeckman et al.l 1)20111) noted that 
such a large momentum injection (ft 10 35 dyne) from a 
small initial radius ro w 100 pc could accelerate a cloud 
with column density Nh ~ 10 21 cm" 2 to a terminal ve- 
locity Voo « 1800 km s -1 . Thus, the energetics of com- 
pact starbursts are sufficient to produce the high- velocity 
outflows we observe, and it is plausible that both radia- 
tion pressure on dust grains and supernova ram pressure 
contribute to driving the winds. 

3.3. Placing These Galaxies in Context 

It is clear from Figure@]that the high-SsFR galaxies in 
our sample constitute a rare population, suggesting that 
they represent an unusual or short-lived phase. Mergers 
of gas-rich galaxies are a viable mechanism for pr oduc- 
ing compact starbursts (jMihos fc Hernquistl [19961) . and 
such gas-rich major mergers are rare at z ~ 0.6 due 
to th e decline in both the gas fraction dTacconi et al 



201C) and the fractional major merger rate ()Lotz et al 
20111) of galaxies since z ~ 2. Furthermore, our high- 



EgpR galaxies are caught in a particular time interval 
where there is both vigorous star formation and strong 
feedback. The length of this phase may be set by the 
gas consumption timescale or the timescale for feedback 
to suppress subsequent star formation. Base d on the 
Kennicutt-Schmidt relation ()Kennicuttl I1998D . a com- 
pact starburst with r e « 100 pc, SFR « 200 M Q , and 
Ssfr ~ 3000 M© yr -1 kpc -2 would have a gas sur- 
face density of S sos ~ 10 11 M© kpc -2 corresponding to 
Mgas ~ 3 x 10 9 Mq inside 100 pc, which would be con- 
sumed on a timescale ~ 20 Myr. This scenario could be 
tested with CO observations of molecular gas masses and 
kinematics for these extreme galaxies. 

Considering our full HST-WISE sample in FigureHJ we 
find values of Ssfr spanning four orders of magnitude. 
This could be explained as an evolutionary sequence 
where the high-SsFR galaxies represent the peak of the 
starburst when the high-velocity outflows are launched, 
while the lower Ssfr galaxies represent a subsequent, 
post-starburst phase. It is interesting to note that all 
12/2 5 galaxies above the S sfr ~ 25 M Q yr" 1 kpc" 2 limit 
from iMeurer et al.1 (|1997h exhibit outflows (with median 
centroid velocity v = —1500 km s" 1 ), while all 7/25 
galaxies without detected outflows (the smallest black 
circles in FigureHJ) have Ssfr < 20 M Q yr" 1 kpc" 2 . If a 



compact starburst is a requirement for the production of 
high-velocity outflows, it may be that the sources with- 
out outflows have not gone through such a phase (see A. 
Robaina, et al., in preparation for a discussion of the re- 
lationship between galaxy morphology, outflow velocity, 
and stellar population age in this sample). 

Finally, we consider the implications of our results for 
models of massive galaxy formation. Simulations of ma- 
jor galaxy mergers with large gas fractions (f ga s ~ 50%) 



can produce M t ~ 10 M Q remnants with 



1 kpc 



(jWuyts et al.ll20l6T ). but our sample includes galaxies of 
similar mass that are smaller in the rest-frame V band 
by almost an order of magnitude. If the mass in these 
galaxies follows their 1^-band light, it would be extremely 
challenging for them to grow in size from r e ~ 0.1 kpc 
to r R ~ 5 kpc to reach the local size-mass relation 



(jShen et al.ll2003h in the t ~ 6 Gyr since z = 0.6. For 
compa rison the compact, quiescent galaxies observ ed at 



npa 

z ~ 2 (jTruiillo et al.ll2007tlvan" Dokkum et al.ll2008t ) have 
t ~ 10 Gyr to grow by a factor of ~ 5. In this context, it 
is worth noting that the half-light radius (ignoring dust 
attenuation) at rest-frame V band can be a factor of 5-10 
smaller than the half-mass radius for a gas-rich merger 
at final coalescence near the peak of starburst activity 
(jWuvts et alj 12010). One could test for such size dis- 
crepancies and probe the radial dependence of the mass- 
to- light ratio for these galaxies by measuring sizes at rest- 
frame near-IR wavelengths. 

4. SUMMARY 

We have measured large SFR surface densities for 
galaxies that exhibit > 1000 km s" 1 outflows. The 
largest Ssfr values are comparable to the Eddington 
limit from radiation pressure on dust grains, and such 
compact starbursts are expected to have substantial mo- 
mentum input from massive stars and supernovae. High- 
velocity outflows have been previously interpreted as a 
signpost of AGN feedback, but given that feedback from 
a compact starburst is capable of producing such a sig- 
nature and is clearly present in this sample, we conclude 
that the outflows we observe are likely driven by star 
formation. 
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